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Context. Radio halos and relics are diffuse radio sources found in galaxy clusters showing significant substructure at X-ray wave- 
lengths . These sources provide important information about non-thermal processes taking place in the intracluster medium (ICM). 
Until now only a few dozen relics and halos are known, while models predict that a much larger number of these sources exist. In this 
paper we present the results of an extensive observing campaign to search for new diffuse radio sources in galaxy clusters. 
Aims. The aim of the observations is to create a large sample of diffuse radio sources in galaxy clusters that help to understand the 
formation of radio relics and halos and can be used to probe the physical conditions of the ICM. 

Methods. We carried out radio continuum observations with the Westerbork Synthese Radio Telescope (WSRT), Giant Metrewave 
Radio Telescope (GMRT) and Very Large Array (VLA) of clusters with diffuse radio emission visible in NVSS and WENSS survey 
images. Optical images were taken with the William Herschel and Isaac Newton Telescope (WHT, INT). 

Results. We discovered 6 new radio relics, including a probable double relic system, and 2 radio halos. In addition, we confirm the 
presence of diffuse radio emission in four galaxy clusters. By constructing a sample of 35 radio relics we find that relics are mostly 
found along the major axis of the X-ray emission from the ICM, while their orientation is perpendicular to this axis. We also compared 
the X-ray luminosity and redshift distributions of clusters with relics to an X-ray selected sample from the NORAS and REFLEX 
surveys. We find tentative evidence for an increase of the cluster's relic fraction with X-ray luminosity and redshift. The major and 
minor axis ratio distribution of the ICM for clusters with relics is broader than that of the NORAS-REFLEX sample. 
Conclusions. The location and orientation of radio relics with respect to the ICM elongation is consistent with the scenario that relics 
trace merger shock waves. 
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1. Introduction 

Radio halos and relics are found in massive merging galaxy clus- 
ters. These radio sources indicate the presence of ma gnetic fields 
and in- situ particle acceleration within the ICM (e.g., Jaffe 1977| 
Govoni & Feretti 2004). Galaxy clusters form through mergers 
with other clusters and galaxy groups as well as through contin- 
uous accretion of gas from the intergalactic medium. Since giant 



dio wavelengths (e.g., |Ensslin et al.||1998HMiniati et~aLl|2000 ). 
The efficiency with which collisionless shocks can accelerate 
particles is unknown and may not be enough to produce the ob- 
served radio brightness of relics. A closely linked scenario is 
that of shock re-acceleration of pre-accelerated electrons in the 
ICM, which is a more efficient mechanism for weak shocks (e.g., 
Markevitch et al.||2005| |Giacintucci et al.|[2008| |Kang & Ryu 
2011b. 



radio halos and relics are found in merging clusters (e.g.,|Buote 



2001 



2009 



Govo ni et al.||2004t |Barrena et al.||2007t |Brunetti et aL 



|Cassano et al.| 2010b), it has been proposed that a small 



fraction of the energy released during a cluster merger event is 
channeled into the (re) acceleration of particles. 

One scenario for the origin of radio relics is that they trace 
merger shock waves within the ICM in which particles are ac- 
celerated by the diffusive shock acceleration (PSA) mechanism 



(Krymskii 1977; Axford et al.| 19771 |Bell| 197 8a b; Blandford&l 
Ostriker 1 97 8) |Druiy| 19 83; Bla ndford & Eichler|1987[lJones~&| 
Ellison|1991[|Malkov & Q'C Drury|2001| ). In the presence of a 
magnetic field these particles emit synchrotron radiation at ra- 



An alternative scenario has been proposed by |Keshet| ( |2010 ) 
based on a secondary cosmic ray electron model, where the time 
evolution of magnetic fields and the cosmic ray distribution are 
taken into account to explain both halos and giant relics. Detailed 
spectral maps, at multiple frequencies, and measurements of the 
magnetic field via the polarization properties can test this model. 
For relics, strong magnetic fields (B > B cm \> = 3.25 x (1 + z) 2 ) 
are predicted and the spectral inde^at the outer edges of relics 
should be flat, with a ^ - 1 . 



F v oc y a , where a is the spectral index 
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Unlike relics, radio halos are found in the center of merging 
galaxy clusters and follow the X-ray emission from the ICM. 
Radio halos have been explained by turbulence injected by re- 
cent merger events. The injected turbulence is thought to be ca- 



pable of re-accelerating relativistic particles (e.g., Brunetti et al. 
|2Q0T|[Petrosian|2001HCassano & Brunetti|2005| ). Alternatively, 
the energetic electrons are secondary products which originate 
from hadronic collisions between relativistic protons and ther- 
mal ions (e.g., |Dennison|[l 980; B lasi & Colafrancesco |"1999; 
Dolag & En Blin 2000; K eshet & Loeb|2010||EnBlin et al. 



2011). 



Recent observations put tension on the secondary models (e.g, 
I Donnert et"aLl|2010a|bt |Jeltema & Profumo||2UTTl |Brown &] 
Rudnick 2011; Bonafede et al. 2011). The existence of ultra- 
steep spectrum radio halos is also claimed not to be in agreement 
with these secondary models (Brunetti et al. 2008 ), but they can 
be explained by the turbulent re-acceleration model. However, 
currently only a few of these ultra-steep spectrum radio halos 
are known so more observations, such as presented in this paper, 
are needed to increase this number and provide better measure- 
ments of the radio spectra. 

In the last decade a number of successful searches have been 
carri ed out to find new diff use radio sources in galaxy clusters 
(e.g. IGiovannini et al.|1999l|Giovannini & Fe retti 2000; Govoml 
etal. 2001 ; Kemp n^&Sarazin|2001HBaccMetal.|2003i|Wn turi 
et al'jpOOT "" 



2008 ; Giovannini et al.||2009[ |van Weeren et~aL 



2009b [ |Rudnick & Lemmerman |2009 ). However, our under 



standing of the formation of these sources is still limited. Models 
for the formation of relics and halos can be tested through statis- 
tical studies of correlations between the non-thermal radio emis- 
sion and global properties of the clusters, such as mass, temper- 
ature, and dynamical state ([Liang et al.|2000| Feretti et al.|2 006 ; 



ICassano et al.|2006||2007||2008H2010a> ." 

We recently discovered two large radio relics in the 



NVSS and WENSS surveys (va n Weeren et al.||20T0| |20TTa> . 
Interestingly, these relics remained unnoticed for about 15 years. 
This suggests that more diffuse radi o sources could be dis- 
covered by ins pection of the NV SS (Con don et al. ||1998[) and 



WENSS (Rengelink et al. 1997) survey images. We therefore 
carried out a visual inspection of the NVSS and WENSS images 
around known cluste rs detec ted by RQSAT ([Voges et al.||1999 



2000 



2007 



Ebeling et al. 



Ebeling et al. 



1998; Bohringer et al. 2000; Kocevski et al. 
2002[|Bohringer et al.|2004| ). 



Candidate clusters hosting diffuse radio emission were ob- 
served with the WSRT, GMRT and/or VLA. Clusters with ex- 
isting published observations were not re-observed. In this paper 
we present the radio images and global properties of the clusters. 
In addition, we investigate the position and orientation of radio 
relics with respect the ICM, and compare the X-ray luminos- 
ity and redshift distributions of clusters with relics to an X-ray 
selected sample. In a follow-up paper we will present the polar- 
ization and detailed spectral properties of the radio emission in 
these clusters. The layout of this paper is as follows. In Sect. [2] 
we give an overview of the observations and the data reduction. 
In Sect. [3] we present the radio images and give an short overview 
of the cluster's properties. We end with discussions and conclu- 
sions in Sects. [4] and \5\ 

Throughout this paper we assume a ACDM cosmology with 
H = 71 km s _1 Mpc -1 , Q m = 0.27, and £l A = 0.73. All images 
are in the J2000 coordinate system. 



2. Observations & Data Reduction 

2.1. Radio observations 

Most of candidate clusters with diffuse radio emission were ob- 
served with the WSRT. GMRT or VLA observations were taken 
of a few clusters which were missed by the WSRT observations. 
A summary of the observations is given in Table [T] The WSRT 
observations were taken in frequency switching mode, alternat- 
ing every 5 minutes between the 21 and 18 cm bands. In this 
paper we will only use the 21 cm data, except for the clusters 
Abell 697 and Abell 2061. The other data will be presented in 
a future paper which will focus on the polarization properties 
and rotation measure synthesis. The total integration time for the 
WSRT observations was 6 hr per cluster, except for Abell 3365 
(see Table [TJ. VLA observations of Abell 523 were taken in 
D-array. We also included archival observations from project 
AB1180 (L-band D and C-array). Abell 3365 L-band VLA ob- 
servations were taken in DnC and CnB array (project AR690). 

GMRT 325 MHz observations (with the hardware correla- 
tor) were taken of Abell 1612 on May 13, 2009, recording both 
RR and LL polarizations with 32 MHz bandwidth. We observed 
Abell 1612 and CIZA J0649.3+1801 at 610/241 MHz with the 
GMRT. The 610/241 MHz observations were taken in dual- 
frequency mode, recording RR polarization at 610 MHz and LL 
polarization at 241 MHz. Total recorded bandwidth was 32 MHz 
at 610 MH z and 6 MHz at 2 41 MHz. The GMRT software back- 
end (GSB; |Roy et al.|2QTQ] ) was used. 

For the data reduction we used the NRAO Astronomical 
Image Processing System (AIPS) package. Standard gain cali- 
bration and RFI removal were performed. Bandpass calibration 
was carried out for observations done in spectral line mode. 
For the 241 MHz observations RFI was fitted and subtracted 
from the data using the technique described by Athreya (2009) 
which was implemented in Obit (Cotton 2008). The fluxes of the 
primary calibrators we re set according to the Perley & Taylor 
(1999) extension to the Baa rs et aT| ( |1977| ) scale. Several rounds 
of phase self-calibration were performed before doing two final 
rounds of amplitude and phase self-calibration. 

For GMRT data in the imaging step we used th e polyhedron 
method ( [Perley |[T989} |Cornwell & Perley | [T992] ) to minimize 
the effects of non-coplanar baselines. Images were made using 
"Briggs" weighting (with robust set to 0.5, Briggs 1995) and 
cleaned down to 2 times the rms noise level (2<x rms ) using clean 
boxes. Finally, we corrected the image for the primary beam at- 
tenuation. For more details on the d ata reduction the reader is 
referred to | van Weeren et aL] ( |201 lb| ). 

2.2. Optical WHT& INT images 

For the clusters Abell 1612, Abell 523, Abell 3365, Abell 2034, 
CIZA J0649.3+1801 and RX J0107.8+5408 we made use of op- 
tical V, R and I band images taken with the WHT (PFIP camera) 
and INT (WFC camera) telescopes between 15-19 April (WHT) 
and 1-8 October (INT), 2009 (for more details see |van Wee ren 
|etal.|2011b| ). 



3. Results 

In this section we present the radio continuum images of the 
clusters. We briefly discuss the results of the radio observations 
for each cluster, a summary is given in Table [2] To compute the 
integrated fluxes for the diffuse radio sources we subtracted the 
flux contribution from the discrete sources. We alphabetically 
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Table 1. Observations 



Cluster 


observation date 


frequency 


bandwidth 


integration time 


map resolution 




rms noise 






MHz 


MHz 


hr 


arcsec 




yt/Jy beam -1 


Abell 1612 


GMRT, 13 May, 2009 


325 


32 


4.0 


11.6" x9.5" 




236 




GMRT, 22 Nov, 2009 


610, 241 


32,6 


2.5,2.5 


7.7" x 4.7", 21" 


xl2" 


64, 777 


Abell 523 


VLA Aug, 25, 2005; 

Nov 28, 2005; Dec 28, 2009 


1425 


37.5 


l.l fll ,3.7 a2 


21" x20" J 




40 


Abell 697 e 


WSRT, 24 Aug, 2009 


1382, 1714 


160, 160 


6.0, 6.0 


34" x2V b 




24,32 


Abell 3365 e 


WSRT, 22 Feb, 2009 


1382 


160 


12.0 


108" x 13" 




29 


Abell 3365 


VLA, 30 Sep, 2009 
30 May, 2009 


1425 


2x50 


0.7 gl ,2.5* 2 


43" x 35", 13.5 


" x 9.2" 


84,27 


Abell 746 e 


WSRT, 19 Sep, 2009 


1382 


160 


6.0 


23" x 18" 




28 


Abell 2034-f 


WSRT, 26 Jul, 2009 


1382 


160 


6.0 


30" x 16" 




24 


Abell 2061' 


WSRT, 23 Jul, 2009 


1382, 1714 


160, 160 


6.0, 6.0 


32" x 16"* 




22,26 


RXC J1053.7+5452 e 


WSRT, 14 Mar, 2009 


1382 


160 


6.0 


24" xl8" 




30 


CIZAJ0649.3+1801 


GMRT, 22 Nov, 2009 


610, 241 


32,6 


3.0, 3.0 


25" x 25" c , 17" 


xl4" 


515, 1800 


RXJ0107.8+5408 e 


WSRT, 29 Aug, 2009 


1382 


160 


6.0 


21" xl7" 




29 



al D array 

al C array 

b the WSRT 1714 MHz image was convolved to the same resolution as the 1382 MHz image 

c convolved to a beam of 25" x 25" 

d D and C array data were combined 

e The minimum baseline length is 36 m. 

f The minimum baseline length is 54 m. 

gl DnC array 

82 CnB array 



labeled these discrete sources for each cluster where applicable. 
The fluxes of the discrete sources are reported in Table [3] and 
they are measured from images made with uniform weighting. 
We included the uncertainties in the subtraction of the discrete 
sources in the uncertainty for the integrated flux measurements 
of the diffuse cluster emission in Tabled We assume that the un- 
certainties for the discrete sources are uncorrelated. More details 
are given in the subsections of the individual clusters. 

We also display overlays onto ROSAT images and galaxy 
iso-density contours (mostly derived from SDSS DR7 photo- 
metric redshifts, Abazajian et al. 2009). The iso-densities were 
computed by counting the number of galaxies per unit surface 
area within a certain photometric redshift range (for SDSS) or 
V-R and R-I color range (for INT images). For the INT im- 
ages we created a catalog of objects using Sextractor (Bertin 
& Arnouts 1996). We then removed all point-like objects (i.e., 
stars) from the catalogs. To exclude galaxies not belonging to 
the cluster we selected only galaxies with R-I and V-R colors 
within 0. 15 magnitude from the average color of a massive ellip- 
tical cD galaxy in the cluster. The range of 0.15 in the colors was 
taken to maximize the contrast of the cluster with respect to the 
fore and background galaxies in the field, but not being too re- 
strictive so that a sufficient number of candidate cluster members 
was selected. 



3.1. Abell 1612 

Abell 1612 is a little studied cluster at z = 0.179 (Popesso 



et al.|2004| ) with a m oderate X-ray lum inosity of L x 01-24 kev = 
2.41 x 1(T 4 erg s" 1 JBohringer et"all|2004| ). In the NVSS sur- 
vey we found an elongated radio source located about 5' to the 
south of the cluster center. The source was completely resol ved 
out in the 1.4 GHz 5" FIRST survey ( [Becker et aL][T995] ) in- 
dicating diffuse emission on scales of about 4'. In our GMRT 



325 MHz image (Fig. [T] left panel) the source has a total ex- 
tent of 4.5', which corresponds to a physical size of 780 kpc. We 
could not identify an optical counterpart for the source (Fig. [3] 
left panel). We therefore classify the source as a radio relic. In 
the 610 MHz image (Fig.[2| left panel) the source is connected 
to a tailed radio galaxy to the north. Combining flux measure- 
ments at 1.4 GHz (NVSS) and GMRT fluxes at 241, 325 and 
610 MHz gives a « -1.4. However, we note that the individual 
flux measurements give a large scatter (Fig.[3j right panel). We 
checked the absolute flux calibration between the different fre- 
quencies by measuring the integrated fluxes of several compact 
sources. This did not reveal any problems. Some of the short 
baselines in the 325 and 241 MHz observations were affected by 
RFI which could have affected the flux measurement since the 
source is quite extended. In addition, the declination is close to 
0° which gives a non-optimal uv-coverage. 

Probably the relic traces a shock in which particles are ac- 
celerated or re-accelerated by a merger-related shock wave. The 
large extent makes it unlikely that the source traces compressed 
fossil radio plasma (Clarke & Ensslin 2006). The connection 
with the tailed radio source to the north favors a re- acceleration 
scenario. In this case, the seed relativistic electrons could be 
supplied by the tailed-radio source. The elongated galaxy iso- 
density contours and ICM, see Fig. [T] (right panel), also hint at a 
merger event along a northwest- southeast axis. 

3.2. Abell 746 

Abel 746 is a Bautz-Morgan (B/M) class III type cluster lo- 

0.01 ([Koester et aT] 



2007). Based 



1 we calculate 



cated at z P hot = 0.232 

on the ROSAT flux from pohringer et al.| ( |2nnfl 

Lx, 0.1-2.4 kev ~ 3.68 x 10 44 erg s" 1 . We discovered an elongated 

radio source located northwest of the cluster center in the NVSS 

survey. The diffuse source is resolved out in the FIRST image, 
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Table 2. Cluster properties 



Relic 


z 


S 1382 MHz 


Pi A GHz 


T d 
Lx, 0.1-2.4 keV 


LLS fl 


classification^ 






mJy 


lO^WHz" 1 


10 44 erg s" 1 


kpc 




Abell 1612 


0.179 


62.8 + 2.6* 


7.9 


2.41 


780 


R 


Abell 746 (center, periphery) 


0.232 


18 + 4,24.5 + 2.0 


3.8,6.8 


3.68 


850,1100 


H, R 


Abell 523 


0.10 


61 +7 / 


1.7 


0.89 


1350 


HorR 


Abell 697 


0.282 


5.2 + 0.5, 4.0 + 0.5* 


2.2 


19.42 


750 


H 


Abell 2061 


0.0784 


27.6 + 1.0, 21.2 ±2.1* 


0.45 


3.95 


675 


R 


Abell 3365 


0.0926 


42.4 + 3.0 J , 5.4 + 0.5* 
42.6 + 2.6^', 5.3 + 0.5^ 


LI'", 0.14* 


0.859 


560^, 235* 


DR 


CIZAJ0649.3+1801 


0.064 


321 + 46< 


\2 h 


2.38 


800 


R 


RX J0107.8+5408 c 


0.1066 


55 + 5 


1.8 


5.42 


1100 


H 


Abell 2034 (center, periphery) 


0.113 


7.3 + 2.0, 24 + 2 


0.28, 0.89, 


3.56 


600, 220 


H or R, R 


RXC J1053.7+5452 


0.0704 


15 + 2 


0.2 


0.44 


600 


R 



a largest linear size 

b R = radio relic, H = radio halo, DR = double radio relic 

c also known as CIZA J0107.7 +5408 or ZwCl 0104.9+ 5350 

d Abell 3365 and Ab ell 1612: Bohri nger et aT] (2004), CIZA J0649.3+1801 and RX J01 07.8+5408 :|Ebeling et al] ( |2002|>, RXC J1053.7+5452 : 

|Popesso etaT] ([2004), Abell 697: ( [Popesso et al.|2004) , Abell 2034 and Abell 2061: |Ebeling et al.| ( [T998| ), Abell 523: |Bohringer eTaLlpOOO] ), 

Abell 746: this work 

e NVSS flux 

/ VLA flux at 1425 MHz 

^ WSRT flux at 1714 MHz 

h using a typical relic spectral index of - 

' GMRT flux at 610 MHz 

j eastern relic 

* western relic 



1.3 



Table 3. Fluxes of compact sources embedded in the diffuse emission 





A 


B 


C 


D 


E 


F 


H I J 


Fi 


gure 


A746 


1.52 + 0.09 


0.36 ± 0.06 


0.24 ± 0.06 


0.37 ± 0.06 


0.45 ± 0.07 


0.93 + 0.06 


0.51+0.07 0.33 + 0.06 0.34 + 0.06 


4 




A523 


10.9 ± 0.2 


5.7 + 0.1 


2.9 + 0.1 










6 




A697 a 


0.50 ± 0.07 


0.29 + 0.06 


0.22 ± 0.06 


0.25 ± 0.08 


0.16 + 0.06 


0.19 + 0.05 




9 




A697^ 


0.38 ± 0.08 


0.24 ± 0.07 


0.17 + 0.05 


0.12±0.05 J 


0.15±0.06 J 


0.07 ± 0.05 J 




9 




A2061 a 


0.53 ± 0.08 


0.41+0.07 


0.36 ± 0.07 










13 




A2061^ 


0.31+0.08^ 


0.33 ± 0.08 


0.20 ± om d 










13 


A3365 c 


1.5 + 0.1 


0.53 ± 0.06 


0.62 ± 0.06 










15 


CIZA J0649.3+ 1801 


7.7 + 0.8 


1.2 + 0.1 


2.5 ± 0.3 










17 


RX J0107.8+5408 


2.6 ± 0.2 


1.1+0.1 


1.0 + 0.1 


0.7 + 0.1 








18 


A2034 


0.93 ± 0.06 


0.59 ± 0.06 


0.34 ± 0.05 










19 


RXC J1053.7+5452 


1.2 + 0.1 


1.4 + 0.1 












2T 



Note: reported fluxes are in mJy 

a \ 382 MHz 

* 1714 MHz 

c from VLA CnB array image 

d by measuring the flux at the 1382 MHz source position 



except for an unresolved source located at the southern tip of the 
diffuse source. In our WSRT image the elongated source is eas- 
ily detected (Fig. [4] left panel) . The source has a largest angular 
extent of 5', which corresponds to a physical size of 1.1 Mpc at 
the distance of Abell 746. The source has a width of 345 kpc and 
is located at a distance of 1 .7 Mpc from the cluster center. 

We do not find it very likely that the point source at the 
southern end of the elongated source is the hotspot of a radio 
galaxy because the second hotspot is missing. In addition, we do 
not detect any radio core. The compact source at the southern 
end of the diffuse source does not have an optical counterpart in 
SDSS images. The SDSS images are contaminated by the bright 
star f Uma (V mag = 4.5). We only conclude that a possible opti- 
cal counterpart should be located farther away than A746. This 



limit on the redshift does not provide any useful constraints to 
rule it out as a giant radio galaxy on the basis of its physical 
size or radio luminosity. The morphology of the source is more 
typical for a radio relic. The size of the sources an its location 
would also agree with this interpretation. As an additional check, 
we analyzed the polarization data from the WSRT observation. 
Radio relics are often polarized at a level of ~ 20% or more (e.g., 
Clar ke & Ensslin|2006[ |Bonafedeetal.|2009b[ |van Weeren etaLl 
2010J). Indeed, we find that the elongated source is polarized up 
to the ~ 50% level (Fig. [5]), which provides additional support 
for the classification as radio relic. 



Galaxy iso-density contours show the cluster to be somewhat 
elongated along a northwest- southeast axis, see Fig. [4] (right 
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Fig. 1. Left: GMRT 325 MHz image of Abell 1612. Radio contours are drawn at levels of [1, 2, 4, . . .] x 4<x rms . Dashed contours are 
drawn at -3<x rms . Right: A1612 X-ray emission from ROSAT, tracing the thermal ICM, is shown by the color image. The original 
image from the ROSAT All Sky Survey was convolved with a 270"FWHM Gaussian. Solid contours are from the GMRT 325 MHz 
image and drawn at levels of [1, 2, 4, . . .] x 3<x rms . Dashed contours show the galaxy iso-density distribution derived from the SDSS 
survey. Contours are drawn at [1.0, 1.4, 1.8, . . .] x 1.1 galaxies arcmin -2 selecting only galaxies with 0.16 < z P hot < 0-20. 
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Fig. 2. GMRT 610 (left) and 241 MHz (right) image of A1612. Contour levels are drawn as in Fig. [l| Black dotted lines in the 
610 MHz image indicate the integration area for the flux measurements. 



panel). The ROSAT image reveals little structure. However, the 
number of photons detected from the cluster is low. 

The WSRT image reveals additional diffuse emission in the 
center of the cluster, albeit at a low SNR. The emission has an ex- 
tent of at least 850 kpc which is typical for a giant radio halo. The 
flux of this radio halo is difficult to estimate due to the low SNR, 
but subtracting the contribution from the compact sources we 
find a flux of about 18 + 4 mJy for the radio halo. To better image 



the diffuse emission we subtracted the clean components from 
the compact sources using an image made with uniform weight- 
ing. We then convolved the image (made with natural weighting) 
to a resolution of Y . The contours from this image are overlaid in 
Fig.|4](right panel). The halo is now better detected and the radio 
emission roughly follows the X-ray emission. The radio power 
of 3.8 x 10 24 W Hz -1 is above the Lx-Pi.4gh z correlation for 
giant radio halos (e.g., Liang et al. 2000[ |Cassano et al. 2006). 



R. J. van Weeren et al.: Radio continuum observations of new radio halos and relics from the NVSS and WENSS surveys 




1.00 








I ' 


; 






*^\ 


i 




- 


0.10 










_ 










P\ 


\j 


001 








i 


- 



1000 
Frequency [MHz] 



Fig. 3. Left: WHT V, R, I color image of Abell 1612. Overlaid are the radio contours from Fig. |2J (left panel). Right: A1612 radio 
relic spectrum. Flux measurements at 241, 325, and 610 MHz are from the GMRT observations. The 1.4 GHz flux is from the NVSS 
survey image. 
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Fig. 4. Left: WSRT 1382 MHz image of A746. Contour levels are drawn as in Fig.[l] Black dotted lines indicate the integration areas 
for the flux measurements. Discrete sources embedded in the diffuse emission are alphabetically labeled, see Table [3] Right: A746 
X-ray emission from ROSAT in orange. The original image from the ROSAT All Sky Survey was convolved with a 225" FWHM 
Gaussian. Solid contours are from the WSRT 1382 MHz natural weighted image convolved to a resolution of 60". Compact sources 
were subtracted and contours are drawn at levels of [1, 2, 4, 8, . . .] x 0.4 mJy beam -1 . Dashed contours show the galaxy iso-density 
distribution derived from the SDSS survey. Contours are drawn at [1.0, 1.4, 1.8, . . .] x 0.6 galaxies arcmin -2 selecting only galaxies 
with 0.16 < Zphot < 0.29. 



The fitted relation from Cassan o et al. ([2006 ) gives a power of 
0.64 x 10 24 W Hz" 1 . We note that both the integrated radio flux 
and L x (from ROSAT) are uncertain, and t hat the intrinsic scat - 
ter in the Lx-^i^ghz relation is quite large (Brunetti et al. 2009). 



The measured radio power is therefore still marginally consistent 
with being on the Lx-.Pi.4GHz correlation. 



3.3. Abell 523 

Abel 523 is a little stu died galaxy cluster located at z = 0.10 
( Struble&Rood|l999|) , with L x , 0.1-2.4 key = 0.89 x 10 44 erg s" 1 
(Bohringer et al. 2000). Galaxy iso-densities derived from INT 
images show a north- south elongated cluster consisting of two 
galaxy clumps. The VLA image (Fig. |6j left panel) reveals a 
large, irregular and diffuse radio source in the cluster as well as 
a number of compact sources related to AGN activity. Our radio 



R. J. van Weeren et al.: Radio continuum observations of new radio halos and relics from the NVSS and WENSS surveys 





1 I ' '0 


1 r ' ' ' 'o ■!? ' ' ' ' 

o o 


40 





- ft' ° * ' 


35 




« o 

5" 

0° 
o o 




cS °-© c ' . '-^^ „_ 


30 


• . o o U — ° ° 


1° 25' 


_0 


^ @ o 








0, i , , ® _ ° - 



00 s 



09'" 30 s 
Right Ascension 



00 s 



Fig. 5. WSRT 1382 MHz polarization map of A746. Total po- 
larized intensity is shown as a grayscale image. Vectors de- 
pict the polarization E- vectors, their length represents the po- 
larization fraction. The length of the E-vectors are corrected 
for Ricean bias (Wardle & Kronberg 1974| ). A reference vec- 
tor for a polarization fraction of 100% is shown in the bottom 
left corner. No vectors were drawn for pixels with a SNR < 4 in 
the total polarized intensity image. Contour levels are drawn at 
[1, 4, 16, 64, . . .] x 40-^ and are from the Stokes I image. 



image and galaxy distribution agree with the recent results from 
Giovannini et al.| (|2011|) which showed the presence of diffuse 



radio emission in this cluster. 

The brightest radio source is a tailed radio galaxy projected 
relatively close to the cluster center. The color of the optical 
counterpart is similar to other galaxies in the cluster making it 
likely that the radio source is associated with cluster. The optical 
INT image is shown in Fig. [7] We also detect some radio emis- 
sion from the largest cD galaxy in the cluster located north of the 
tailed radio source. In the southern part of the cluster there is a 
brighter compact source associated with another large elliptical 
galaxy. 

The diffuse source has a patchy morphology, with the bright- 
est part of the diffuse source located to the northwest of the tailed 
radio source. To the west the diffuse source extends into two fil- 
amentary structures. The total flux of the diffuse radio source, 
minus the point sources and head-tail galaxy, is 61 ± 7 mJy. The 
diffuse source has a largest extent of 1.35 Mpc. Both numbers 
are consistent with the result from Giovannini et al. ('201 1). 

The large extent and morphology make it unlikely that the 
diffuse sources is direct ly related to the tailed radio galaxy. 
|Giovannini et al.| ( |2011| ) classify the source as a radio halo. 
They note that the radio emission permeates both galaxy clumps. 
However, the possibility that the source is a radio relic projected 
onto central region of the cluster should also be considered. The 
source has a very patchy morphology unlike typical radio halos 
for which the surface brightness follows that of the X-ray emis- 
sion. Currently we have too little constraints to completely rule 
out the relic scenario. Deep X-ray observations will be needed to 
characterize the dynamical state of the cluster. Polarization ob- 
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Fig. 7. INT V, R, I color image of Abell 523. Overlaid are the 
radio contours from Fig. [6] 



servations can be used to distinguish between the radio halo and 
relic scenarios as radio relics are usually strongly polarized (at 
the 10-20% level or more) while halos are mostly unpolarized. 

3.4. Abell 697 

Abell 697 is a massive Bautz-Morgan (B/M) type II-III clus- 
ter, located at z = 0.282, with a hig h X-ray luminosity o f 



L x , 0.1-2.4 kev = 19.42 x 10 44 erg s _1 (Popesso et al. 2004) 



The ICM has an elliptical shape and an overall temperatu re of 
8.8^ keV determined from Chandra observations by Maughan 
|et al.|p008t . |Cavagnolo et al.|p009|) reported a slightly higher 
temperature of 9.52 keV. |Girardi et al. ( |2006| ) found a velocity 
dispersion of 1334^^ 4 km s" 1 for the cluster. They noted this 
is expected in the case of energy-density equipartition between 
galaxies and gas. They suggested that the cluster has undergone a 
complex cluster merger event, occurring mainly along the LOS, 
with a transverse component in the SSE-NNW direction. 

|Kempner & Sarazrn] ( |2001| ) first suggested the presence of 
diffuse radio emission in the cluster. [Venturi et al. (2008 ) showed 
the presence of a radio halo usi ng 610 M Hz GMRT o bservations, 
the halo was also reported b y |Rudnick & Lemmerman| p009| ) 
and [Giovannini et aT] ([2009 ). Ma cario et al. ( [201 0| ) presented a 
more detailed study which included GMRT 325 MHz observa- 
tions. They found the radio halo to have a an ultra- steep spectrum 
(^32 4 5 ™Hz) with a about -1.7 to -1.8. 

We detect the radio halo at both 1382 and 1714 MHz with 
better SNR than the previous observations at 1.4 GHz. In our 
1382 MHz image we find a total extent of about 750 kpc fo r the 
radio halo w hich is lower than the 1.3 Mpc reported by Macario 
|et~aLl ( |20T0l ) at 325 MHz. The lower 1382 MHz extent is ex- 
pected due to the steep radio spectrum of the halo. To better 
image the diffuse emission and remove the contribution from 
compact sources, we made images with uniform weighting and 
excluded data < 2.5kA at both 1382 and 1714 MHz. We then sub- 
tracted the clean components of the compact sources from the 
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Fig. 6. Left: A523 VLA 1.4 GHz image. Contour levels are drawn as in Fig. [l| Black dotted lines indicate the integration area 
for the flux measurement. Discrete sources embedded in the diffuse emission are alphabetically labeled, see Table [3] Right: A523 
X-ray emission from in orange. The original image from the ROSAT All Sky Survey was convolved with a 225" FWHM Gaussian. 
Solid contours are from the VLA 1425 MHz image and drawn at levels of [1,2, 4, 8, . . .] x 4<x rms . Dashed contours show the galaxy 
iso-density distribution derived from INT images. Contours are drawn at [1.0, 1.2, 1.4, . . .] x 1.6 galaxies arcmin" 2 selecting only 
galaxies with colors 1.12 < V - R < 1.42, 0.25 < R - I < 0.55, i.e., within 0.15 magnitudes the V-R and R-I color of the central 
cD galaxy. 



uv-data before re-imaging. To increase the SNR we combined 
the 21 and 18 cm images after convolving them to a common 
resolution of 29" x 17", see Fig.[9] 

The discrete sources D, E, and F (see Fig. [9]) were not 
detected with a SNR > 3. By using the positions from the 
1382 MHz image, we could still get an approximate estimate 
of their fluxes. We subtracted the flux fro m the discrete s ources 
(Table [3]) and took the fluxes reported by Maca rio et al. ( 2010| ) 
to fit a power-law radio spectrum through the flux measurement 
with a = -1.64 ± 0.06, see Fig. [8 This confirms that the radio 
halo has a very steep spectral index, a lthough it is margina lly 
flatter than the -1.7 to -1.8 reported by |Macario et al.|p010) . 

We do not detect any polarized flux from the radio halo in 
our WSRT observations. We set an upper limit on the radio halo 
polarization fraction of 6% at 1382 MHz. 



3.5. Abell2061 

Abell 2061 is a Bautz-Morgan (B/M) type III cluster lo- 
cated at z = 0.0784. The cluster has a X-ray lum inosity of 
L x , 0.1-2.4 kev = 3.95 x 10 44 erg s" 1 (Eb eling et all 1998). The 



cluster h as a velocity dispersion of 673^ km s" 1 dOegerle & 
Hill|2001 ). The ROSAT image shows the cluster to be extended 
in the northeast- southwest direction, see Fig.[l0j Previous stud- 
ies report that the cluster is part of the Corona Borealis su- 
percluster ( [Small et aL]|1998[ |Marini et a"T]|2004| ). The galaxy 
iso-density contours display three substructures located roughly 
along a northeast southwest axis. 

Abell 2067 is located at a projected distance of only 
1.8 h~ l Mpc north from Abell 2061, and the clust er's system- 
atic v elocities are separated by ~ 1600 km s _1 ( Qegerle & Hill 
2001| ). Therefore, they probably form a bound system ( Marini 



et al.||2004| |Rines & Diaferio||2006| ) consisting of a massive 
cluster (A2061) with a smaller in-falling group/cluster (A2067). 
A2061 also contains an X-ray extension in the direction of 
A2067 (towards the north-east) which also sugges ts a dynam - 
ical connection between the two systems (Ma rini et al.||2004] ). 
According to Mari ni et"aL} the interaction between the clusters 
is in the phase in which the cores have not yet met and in which 
the formation of a shock is expected. The global temperature for 
A2061 is reported to be 4.53^- 4 g keV from BeppoSax observa- 
tions (Marini et al. 2004). A region with a higher temperature 
is found in the northern part of A2061 with a temperature of 
10.67+247 keV' This region could correspond to the presence of 
an internal shock. 



Kempner & Sarazin| ( |2001| ) reported a possible relic in the 



southwest periphery of A2061 in WENSS images. They found 
a flux of 104 ± 15 and 19 + 3 mJy at 327 MHz and 1.4 GHz, 
respectively. This would give a spectral index of a = -1.17 + 
0.23. Rudnick & Lemmerman (2009) also listed the presence of 
the diffuse peripheral source and measured a flux of 120 mJy in 
the WENSS image. They also claimed the presence of additional 
diffuse emission in the center of the cluster which could make up 
a radio halo. 

With our WSRT observations we confirm the presence of the 
diffuse radio source in the southwestern periphery of A2061 and 
we also classify it as radio relic. We cannot confirm the presence 
of a radio halo. The radio relic is clearly seen in both the 21 and 
18 cm WSRT images, Fig. 



11 



The relic is located at a distance 
of 1.5 Mpc from the cluster center and has a largest angular 
extent of 7.7', corresponding to a physical size of 675 kpc. In 
the direction towards the cluster center the relic has an extent of 
~ 320 kpc. The western outer boundary of the relic is somewhat 
more pronounced, while the emission fades more slowly inwards 
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Fig. 8. Left: WSRT 1382 MHz image of Abel 697. Contour levels are drawn as in Fig.[l| Black dotted lines indicate the integration 
area for the flux measurements. Middle: WSRT 1714 MHz image of Abell 697. Contour levels are drawn as in Fig. [I] Right: A697 
radio halo spectrum. Flux measurements at 325 and 610 MHz are taken from|Macario et al.](|2010|). 
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red color. For comparison, we have marked two galaxies located 
approximately at the distance of A2061 (Fig. 13 ). 

The radio spectrum is fitted by a single power-law spectrum, 
see Fig. [12] We find a = -1.03 + 0.09 . The 325 MHz WENSS 
flux is from Rudnic k & Lemmerm an The 325 MHz measure- 
ment is however uncertain as the SNR on the relic is low. To 
better constrain the spectral index an additional low-frequency 
flux measurement is needed. 

The picture that emerges from our observations is that of a 
radio relic tracing a shock wave from a cluster merger event. The 
merger event does not seem to be directly related to the cluster 
A2067, located north of A2061. Instead the shock wave is more 
likely related to the observed substructures seen in the SDSS 
galaxy iso-density contours. 

3.6. Abell 3365 
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Fig. 9. WSRT combined 1382 and 1714 MHz image of the ra- 
dio halo in Abel 697, compact sources were subtracted from 
the uv-data. Red contour levels are drawn at [1,2,4,8...] x 
60 yuJy beam -1 . Blue contours display the compact sources at 
1382 MHz. This high-resolution image (19" x 10") was made 
with uniform weighti ng and data < 2.5kA was excluded. Contour 
levels are drawn at VU>2, 4, 8. . .] x 100 //Jy beam -1 . Discrete 
sources embedded in the diffuse emission are alphabetically la- 
beled, see Table [3] 



to the cluster center. The relic consist of a northern and fainter 
southern component. Two compact sources are found directly to 
the north and south of the relic. 

An optical color image at the location of the relic is shown 
in Fig. 13 This image does not reveal any obvious counterparts 
to the radio relic. One galaxy is located at the brightest region 
of the relic, but no compact radio source is associated with that 
galaxy. The bright compact sources to the north and south of 
the relic are both associated with background galaxies unrelated 
to A2061 because of their small angular sizes (< 3") and very 



Abell 3 365 ([Abell et al |[l989) is located at a redshift of z = 
0.0926 flStruble & Rood|1999| ). The cluster is listed as a Bautz- 
Morgan (B/M) t ype II cluster. The E SQ Nearby Abell Cluster 
Survey (ENACS, Katgert et al. 1996) measured a velocity dis- 
persion of 1153 km s" 1 ( Mazure et al.||l996 ) for A3 3 65, quite 
high compared to other clusters in the sample. Galaxy cluster 
RXC J0548.8-2154 is located 9' to the west of the NED listed 
position for A3365 at a redshift of z = 0.0928. The galaxy dis- 
tribution, see Fig. [14] (right panel), peaks at the center of A3365, 
while the X-ray peak is located at the position of RXC J0548.8- 
2154. The galaxy distribution around A3365 is complex, with 
two main concentrations along an east-west axis and a smaller 
concentration at the far west. The ICM distribution is also com- 
plex and elongated in the east- west direction. Based on the very 
similar redshift, and complex X-ray emission and galaxy distri- 
bution, we conclude that A3365 and RXC J0548.8-2154 belong 
to the same merging system to which we will simply refer as 
A3365. 

In the NVSS image we noted the presence of a peripheral 
elongated radio source to the east of the cluster center. The 
WSRT image (Fig. 14 left panel) reveals a second smaller dif- 
fuse source on the west side of the cluster. The diffuse sources 



are also detected in the VLA images (Fig. [T5]). For both sources 
we cannot identify obvious optical counterparts in our INT im- 
ages (Fig. [16]). The east and western sources have an angular 
extent of 5.5' and 2.3', which correspond to 560 and 235 kpc at 
the distance of A3 365. We classify the eastern source as a ra- 
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Fig. 11. Left: A2061 WSRT 1382 MHz image. Contour levels are drawn as in Fig.[l| Black dotted lines indicate the integration area 
for the flux measurements. Right: A2061 WSRT 1714 MHz image. Contour levels are drawn as in Fig.[T] 



50 



45 



40 



35 



30 



30 u 25' -0 









v o \ x a x n 



• 



-ff « \o\ |° fl ^\ 



15 h 22 m 00 s 21 m 30 s 



^ - -U 6- 



$ 

■ °. ■ I ■ 



00 13 
Right Ascension 



20 m 30 s 




1.00 



0.10 




1000 
Frequency [MHz] 



Fig. 12. A2061 radio relic spectrum. Flux measurement at 
325 MHz is taken from Rudnick & Lemmerman (2009). 



Fig. 10. A2061 X-ray emission from ROSAT in orange. The 
original image from the ROSAT All Sky Survey was convolved 
with a 225 " FWHM Gaussian. Solid contours are from the 
WSRT 1382 MHz image and drawn at levels of [1,2, 4, 8, . . .] x 
4<r rms . Dashed contours show the galaxy iso-density distribu- 
tion derived from the SDSS survey. Contours are drawn at 
[1.0, 1.4, 1.8, . . .] x 0.55 galaxies arcmin -2 selecting only galax- 
ies with 0.05 < Zphot < 0.1. 



distribution which suggests a merger event along an east-west 
axis (Fig. [14] right panel). We also made an image where we 
subtracted the compact sources from the uv-data. To make this 
image, we subtracted the clean components from a uniformly 
weighted image (CnB array), from the DnC array data. This im- 
age is overlaid in Fig. [14] 

3.7. CIZAJ0649.3+1801 



dio relic. Most likely, the western source is another radio relic, CIZA J0649.3+1801 is a little studied galaxy cluster located at 

which means that A3365 hosts a double radio relic system, but z = 0.064 discovered by Ebeling et al. (2002) at a galactic lati- 

more observations are needed to confirm this classification. This tude of b = 7.668°. The cluster has a moderate X-ray luminosity 

interpretation is supported by the elongated X-ray and galaxy of Lx, 0.1-2.4 kev = 2.38 x 10 44 erg s" 1 . The cluster forms part of a 
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Fig. 14. Left: A3365 WSRT 1382 MHz image. Contour levels are drawn as in Fig.[j] Right: A3365 X-ray emission from ROSAT 
in orange. The original image from the ROSAT All Sky Survey was convolved with a 225" FWHM Gaussian. Solid contours 
are from a VLA 1.4 GHz image with compact sources subtracted from the uv-data. This DnC array image has a resolution of 
47" x 42". Contours are drawn at levels of [1, 2, 4, 8, . . .]x0.4 mJy beam -1 . Dashed contours show the galaxy iso-density distribution 
derived from INT images. Contours are drawn at [1.0, 1.1, 1.2, . . .] x 0.78 galaxies arcmin -2 selecting only galaxies with colors 
0.6 < V - R < 0.9, 0.54 < R - I < 0.84, i.e., within 0.15 magnitudes the V-R and R-I color of the largest cD galaxy. 
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Fig. 15. A3365 VLA 1.4 GHz DnC (left) and CnB (right) array images. Contour levels are drawn as in Fig. []J Black dotted lines 
indicate the integration areas for the flux measurements. Discrete sources embedded in the diffuse emission are alphabetically 
labeled, see Table [3] 



supercluster in the Zone of Avoidance hidden by the Milky Way 
flKocevski et al.|2007) . 

We discovered the presence of a diffuse ~ 10' elongated 
source to the west of the cluster in the NVSS survey. We do not 
find an optical counterpart for the source in our WHT images. In 
the GMRT image, see Fig. [17] (left panel), the source has a total 
extent of 10.6' which corresponds to a size of about 800 kpc at 



the distance of CIZA J0649.3+1801. We therefore classify the 
source as a peripheral radio relic located at 0.8 Mpc from the 
cluster center. The relic is also visible in the 241 MHz image, 
although the SNR on the relic is very low (Fig. [17] right panel). 



The ROSAT image reveals another fainter X-ray source lo- 
cated between the radio relic and the cluster center. The source 
is not resolved and therefore it could be unrelated to the cluster. 
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Fig. 16. A3365 optical INT V, R, I color images at the location of the radio relics. Radio contours are from Fig. 15 (right panel). 




3.8. CIZAJ01 07.7+5408 
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Fig. 13. A2061 optical WHT V, R, I color image at the location 
of the radio relic. Radio contours are from Fig.[TT] The two red 
circles indicate g alaxies with spectrosc opic redshifts of 0.0815 
(west) and 0.061 (Abazajian et al. 2009). Discrete sources em- 
bedded in the diffuse emission are alphabetically labeled, see 
Tabled 



Also we do not detect any group or cluster of galaxies associ- 
ated with this source in our optical WHT images. Therefore, this 
X-ray source seems to be unrelated to CIZA J0649.3+1801. 



CIZA J0107.7+5408 is located at z = 0.1066 ( [Crawford et aL] 
|1995[ ) and has quit e a high X-ra y luminosity of L x , 0.1-2.4 kev = 
5.42 x 10 44 erg s" 1 JEbeling et al.|2002| ). The cluster is projected 
relatively close to the galactic plane with b = -8.65°. Both the 
NVSS and WENSS survey images display an extended diffuse 
radio source located roughly at the cluster center. Our WSRT im- 
age clearly reveals the presence of a somewhat elongated radio 
halo with a largest extent in the north- south direction of 1 . 1 Mpc. 
The galaxy and ICM distribution are also elongated along the 
major axis of the radio halo, see Fig. 18 which supports the sce- 
nario that the cluster is currently undergoing a merger event. An 
image of the radio halo with the discrete so urce s subtracted (us- 
ing the same technique as described in Sect.|3.4[) is overlaid with 



contours. The radio power of 3.8 x 10 24 W Hz 1 is consistent 
with the Lx-^P 1.4GHz correlation for giant radio halos. 



3.9. Abel 1 2034 

Abell 2034 is a merging galaxy cluster a s revealed by a Chandra 
X-ray study from Kemp ner etaL] ( |2003| ). |White| ( |2000| ) reported 
a global cluster temperature o f 7.6 keV in agreement with the 
value of 7.9 



0.4 keV found by Ke mpner etaL] ([2003 ). Kempner 
& Saraziii] ( |2001| ) suggested the presence of a radio relic on the 



basis of a WENSS image. The diffuse emissi on is located clos e 
to the position of a cold front found by Kempner et al. ( 2003 ). 
The presence of diffuse emission was confirmed by |Giovan nini 



et al.|(2009) and was classified as an irregular elongated radio 



halo. 

We also observe the diffuse emission to brighten at the po- 
sition of the cold front in the WSRT image. The SDSS galaxy 
distribution (Fig. 19 right panel) is bimodal which supports 
the scenario that the cluster is undergoing a merger event. We 
consider the classification of the source uncertain, since the 
radio emission does not show a very clear correlation with 
the X-rays. It could therefore also be a radio relic. The ra- 
dio power CPi.4ghz = 0-28 x 10 24 W Hz" 1 ) though is in agree- 
ment with the Lx-.Pi.4GHz correlation for giant radio halos (e.g., 



R. J. van Weeren et al.: Radio continuum observations of new radio halos and relics from the NVSS and WENSS surveys 



13 



« 18 u 00' 



17 u 50' 




06" 50'" 00 



49 30 3 00 3 

Right Ascension 



' ' i ' / i i i i i i v - i 

(J , y r 

$ £/' - 

'. & 

- 

o 



06" 50'" 15 s OCT 



Right Ascension 



Fig. 17. Left: CIZA J0649.3+1801 X-ray emission from ROSAT in orange. The original image from the ROSAT All Sky Survey 
was convolved with a 225" FWHM Gaussian. Solid contours are from the GMRT 610 MHz image convolved to a resolution of 
25" and drawn at lev els of [1, 2, 4, 8, . . .] x 4cr rms . The source NVSS J064928+ 175700 was removed using the "peeling"-method 
(e.g., Noordam 2004). Black dotted lines indicate the integration area for the flux measurement. Discrete sources embedded in the 
diffuse emission are alphabetically labeled, see Table [31 Right: CIZA J0649.3+1801 GMRT 241 MHz image. Contour levels are 
drawn as in Fig.[l] The dashed line is the 15 mJy bearn^ contour of the 241 MHz images convolved to a circular beam of 45". 
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Fig. 18. Left: CIZA J0107.7+5408 WSRT 1382 MHz image. Contour levels are drawn as in Fig.[j] Black dotted lines indicate 
the integration area for the flux measurement. Discrete sources embedded in the diffuse emission are alphabetically labeled, see 
Table[3] Right: CIZA J0107.7+5408 X-ray emission from ROSAT in orange. The original image from the ROSAT All Sky Survey 
was convolved with a 225" FWHM Gaussian. Solid contours are from the WSRT 1382 MHz image convolved to a resolution of 40". 
Compact sources were subtracted and contours are drawn at levels of [1, 2, 4, 8, . . .] x 0.3 mJy beam -1 . Dashed contours show the 
galaxy iso-density distribution derived from INT images. Contours are drawn at [1.0, 1.2, 1.4, . . .] x 0.3 galaxies arcmin -2 selecting 
only galaxies with colors 0.85 < V - R < 1.15, 0.6 < R - I < 0.9, i.e., within 0.15 magnitudes the V-R and R-I color of the central 
cD galaxy. 
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Fig. 20. A2034 WHT V, R, I color image around the small radio 
relic. Contours at 610 MHz from the GMRT are overlaid and 
drawn at levels of [1, 2, 4, 8, . . .] x 4<x rms . The noise level in the 
image is 41// Jy beam -1 and the beam size is 6.8" x 4.3". 



Liang et aL|2000| |Cassano et aL][20Q6| ). The total flux we find 
is 7.3 + 2.0 mJy, lowe r than the 13.6 ± 1.0 mJy reported by 
|Giovannini et al.|(|2009|). The reas on for the difference is unclear 
because |Giovannini et al. ( |2009 ) do not report the integration 
area and the fluxes of the discrete sources that were subtracted. 
Deeper observations which sufficient resolution are needed to 
classify the nature of the diffuse emission. Polarimetric observa- 
tions could be very helpful here as radio relics are usually highly 
polarized. 

We report the detection of a new relatively small radio relic 
located west of the cluster center (NVSS J150940+333119, 
WN B 1507.6+3342, 7C 150739.39+334252.00). This relic is al- 
ready visible in the previous studies mentioned above but was 
not recognized as a radio relic probably because it is quite com- 
pact. The source has a size about 220 by 75 kpc and a spec- 
tral index « -1.2, including flux measurements from the NVSS, 
WENSS and 7C ( [Hales et al.|2007| ) surveys . A high-resolution 
610 MHz GMRT of the relic overlaid on an optical WHT image 
is shown in Fig. 20 The source does not have an optical counter- 
part which should have easily been visible in the WHT image. 
The radio plasma could have originated from the compact ra- 
dio source located immediately south of the relic. On the other 
hand the spectral index is more typical of a relic directly tracing 
a shock wave. 



3.10. RXCJ1 053.7+5452 



RXC J1053.7+5452 is located at z = 0.0704. Aguerri et al. 



(2007) reported a velocity dispersion of 665^ km s A and 



an r 2 oo radius of 1.52 Mpc based on SDSS data. Rudnick & 
Lemmerman (2009 ) mentioned the presence of a diffuse radio 



reli c with an extent of about 1 Mpc on the west side of the clus- 
Rudnick & Lemmerman| reported a total flux of 0.36 Jy at 



ter. 



325 MHz. The low-surface brightness radio relic is also vis- 
ible in the WSRT image. We find a total extent of 600 kpc, 



lower than that of Rudni cTT&Lemmermanl ( [2009] ). The ROSAT 
image overlaid with SDSS galaxy iso-density contours and the 
1382 MHZ WSRT image are shown in Fig. [21] The galaxy distri- 
bution seems irregular with the main peak to the southeast of the 
X-ray center and a second peak to the northwest. The radio relic 
is roughly located along a line connecting these two galaxy con- 
centrations. This is expected for a cluster merger event, with the 
shock waves traveling outwards along the merger axis. Deeper 
radio observations are necessary to better map the extent of this 
faint radio relic. 

4. Discussion 

Giant radio relics are proposed to directly trace merger shock 
waves in galaxy clusters. Simulations show that in the case of 
a binary merger event two shell-like shock waves form at the 
moment of core passage (e.g., Roettiger et al.|[T99 9; Ric ker & 
Sarazin 2001 ). These shock waves then move outwards into the 
lower density ICM of the cluster's outskirts. Double radio relics 
are thought to trace these binary merger events, with two sym- 
metric shock waves on each side of the cluster center. However, 
often merger events are more complex with multiple substruc- 
tures merging, so that relics are not necessarily symmetric struc- 
tures and not always come in pairs. Also, the shock structures 
may break up when they interact with the galaxy filaments con- 
nected to the cluster (e.g., Paul et al. 2010). Still, it is expected 
that relics are mainly found along the main axis of a merger 
event, while their orientation is perpendicular to this axis. 

We have tested this prediction by constructing a sample of 35 
relics, taking relics from the literature and those presented in this 
paper (see TableQ. We did not include any radio relics classified 
as AGN relics or radio phoenices and selected only relics in clus- 
ters which are detected in the ROSAT All Sky Survey (RASS) 
images. 

For every relic we record the end positions of their largest 
spatial extent (Rl, R2), see Fig. 22 The line connecting these 
two points we define as the relic's major axis. The midpoint be- 
tween the two extrema we take as a the relic's center position. 
The RASS image can be used to estimate the position of a pos- 
sible merger axis. We convolve the ROSAT X-ray images for 
each cluster to the same spatial resolution of 650 kpc. We then 
fit a 2-dimensional elliptical Gaussian to the X-ray emission. We 
record both the position of the major axis and the center of the 
fitted Gaussian. For the merger axis we take as a proxy the ma- 
jor axis of the fitted X-ray emission. Finally, we compute the 
angle (a Y ) between the major axis of the ICM and the line clus- 
ter center-radio relic center. The resulting histogram is shown in 
Fig. [22] From this histogram we see that relics are preferably 
found along the major axis of the ICM. This is in line with the 
simple picture that shock waves propagate outwards along the 
merger axis. We also calculate the angle (fl) between the major 
axis of the ICM and the relic's major axis. We find that most 
relics are oriented perpendicular to the ICM major axis, also in 
agreement with a proposed shock origin for radio relics. 



4.1. Comparison with the REFLEX and NORAS X-ray 
clusters 

We also compared the properties of cluster hosting radio relics 
with the X-ray clusters from the NOR AS ( [Bohringer et al.|2000| ) 
and REFLEX ( Bohringer et al.|2004| ) surveys. The NORAS sur- 
vey contains 378 galaxy clusters and has an estimated com- 
pleteness of about 50% at an X-ray Flux of F x , 0.1-2.4 kev = 
3.0 x 10" 12 erg s" 1 cm -2 . The REFLEX sample contains 447 



33 u 25' 



R. J. van Weeren et al.: Radio continuum observations of new radio halos and relics from the NVSS and WENSS surveys 

40 P ' ' 



15 




33" 25' 



30 s 15 s OCT 

Right Ascension 




30 b 15 3 00 3 

Right Ascension 



09'" 45* 



Fig. 19. Left: A2034 WSRT 1382 MHz image. Contour levels are drawn as in Fig.[l| Black dotted lines indicate the integration areas 
for the flux measurements. Discrete sources embedded in the diffuse emission are alphabetically labeled, see Table [3] Right: A2034 
X-ray emission from ROSAT in orange. The original image from the ROSAT All Sky Survey was convolved with a 180" FWHM 
Gaussian. Solid contours are from the WSRT 1382 MHz image and drawn at levels of [1, 2, 4, 8, . . .] x 4<x rms . Dashed contours show 
the galaxy iso-density distribution derived from the SDSS survey. Contours are drawn at [1.0, 1.2, 1.4, . . .] x 1.1 galaxies arcmin -2 
selecting only galaxies with 0.07 < z p hot < 0-16. 



galaxy clusters and has and a completeness of about 90% at the 
same flux level as the NORAS survey. For each cluster in this 
sample, we fitted a 2-dimensional elliptical Gaussian to the X- 
ray emission from the RASS image, using the same procedure as 
for the clusters hosting giant relics. For some clusters the fitting 
procedure did not converge because of nearby bright confusing 
X-ray sources, these cluster were not included in the analysis. 
The resulting histogram is displayed in Fig. [23] The distribu- 
tion of the major-minor axis ratios for the relic cluster sample is 
broader. This is expected since relics should be found in merging 
cluster which are typically more elongated. 

In addition, we compare both the X-ray luminosity and 
redshift distribution of the NORAS-REFLEX sample with the 
relic sample. We selected all clusters with a flux larger than 
3.0 x 10 -12 erg s _1 cm" 2 in the ROSAT band. Below this flux 



limit the NORAS is more than 50% incomplete, see Fig. 24 



The total number of clusters above this flux limit is 540, and the 
number of clusters with relics is 16 (using the same flux cutoff). 
From this we find that the currently observed fraction of clusters 
hosting relics is 3% in this sample. The list of relics is given in 
Table |5] and the radio power are plotted as function of redshift 



are plotted in Fig. 25 



The resulting histograms are displayed in Fig. 26 Although 
the number of known relics is rather small, and the relics 
were selected using various methods, there are few interesting 
trends visible. Apparently, the fraction of clusters with relics 
increases with the X-ray luminosity, from about a percent at 
Lx, 0.1-2.4 kev = 1 x 10 44 erg s" 1 to more than 10% above 
~ 1 x 10 45 erg s" 1 . Also, the redshift distribution for cluster with 
relics is somewhat broader than the corresponding NORAS- 
REFLEX sample. Therefore, the chance of finding a relic above 
a certain flux density increases with redshift for clusters se- 
lected from flux-limited X-ray surveys. Since the number of 



clusters with relics is small this might be a statistical fluctua- 
tion. However, given that the average fraction of clusters with 
relics is 3 %, the probability to find four (or even more) in the 26 
clusters with Lx > 10 45 erg s" 1 is 0.7%. Hence, a pure statistical 
fluctuation is rather unlikely. This rises the question if selection 
effects may causes this trend. 

As discussed above the NORAS-REFLEX sample is more or 
less complete up to 3 x 10" 12 erg s" 1 cm -2 . The relic sample in 
contrast is probably not complete up to a specific flux limit. For 
example, large low-surface brightness relics could be missed, as 
are relics in more distant clusters because they are barely re- 
solved in the NVSS and WENSS surveys images and therefore 
not easily recognizable. The construction of a flux-limited relic 
sample is therefore challenging as selection effects due to angu- 
lar size, morphology and surface brightness need to be properly 
taken into account. In addition, one needs to properly identify 
other diffuse radio sources (such as radio halos and giant radio 
gala xies) to prevent th em from ending up in the relic sample. 
Rudnick & Lemmerman (2009 ) give a good overview of some 



of the problems involved in constructing these samples. The im- 
portant question is whether there is any systematic effect which 
leads to a preferential detection of relics in luminous clusters. 
Since the luminous clusters are on average at higher redshift, 
see Fig. |24j relic detection is less affected by a too low surface 
brightness. However, it is possibly affected by resolution effects. 
This will only decrease the detection probability at higher red- 
shift (luminosity) and hence there is no evident reason why relics 
in more luminous clusters should have a better chance of being 
discovered. 

Skillman et al.| ( [20TT] ) found in simulations, based on DSA, 
that the radio power of clusters with relics scatters largely for a 
given X-ray luminosity. In addition, they found that the mean ra- 
dio power strongly correlates with the X-ray luminosity. Hence, 
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Fig. 21. Left: RXC J1053.7+5452 WSRT 1382 MHz image. Contour levels are drawn as in Fig.Q Black dotted lines indicate the 
integration area for the flux measurement. Right: RXC J1053. 7+5452 X-ray emission from ROSAT in orange. The original image 
from the ROSAT All Sky Survey was convolved with a 225"FWHM Gaussian. Solid contours are from the WSRT 1382 MHz 
image and drawn at levels of [1, 2, 4, 8, . . .] x 4<x rms . Dashed contours show the galaxy iso-density distribution derived from the 
SDSS survey. Contours are drawn at [1.0, 1.2, 1.4, . . .] x 0.25 galaxies arcmin -2 selecting only galaxies with 0.05 < z p hot < 0.1. 
Discrete sources embedded in the diffuse emission are alphabetically labeled, see Table [3] 
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Fig. 22. Left: Schematic illustration of the angle between the major axis of the ICM and the line relic center-cluster center (a T ), 
and the angle between the relic orientation and major axis of the ICM (J3), see Table [4] Top Right: Histogram of angles between the 
major axis of the X-ray emission and the line connecting the cluster center with the center of the relic. Bottom Right: Histogram of 
angles between the major axis of the X-ray emission and the relic's major axis. 



the fraction of clusters with relics should increase with X-ray lu- 
minosity. To compute the actual fractions many factors need to 
be considered, namely the X-ray flux limit, the discovery prob- 
ability for radio relics, the radio power distribution of relics and 
the abundance of clusters as function of X-ray luminosity and 
redshift. In Nuz aet al. (201 1 ) we cary out this analysis in detail. 
We postulate a radio relic probability density based on numeri- 
cal simulations p(P, M v i r , z, v 0D sX where P is the radio luminos- 
ity, M v i r is the virial mass of the cluster and v 0D s is the observing 
frequency. We convolve this with the cosmological abundance 



of dark matter halos. As a result we indeed find that the fraction 
of radio relics in the NORAS -REFLEX sample should increase 
with both the X-ray luminosity and the redshift, see Fig. 26 
(solid lines). 

The reason why in simulations more massive clusters show 
on average much brighter radio relics is unclear. We speculate 
that multiple aspects contribute: temperature and density are 
higher and shock fronts are larger in more massive clusters. It 
seems also likely that mergers of more massive clusters result in 
higher Mach numbers. Since according to the sub-grid models 
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Cluster name 



a r 



"location PA" (deg) "orientation PA" (deg) 



Abellll5 0.197 

Abell 521 0.247 

Abell 523° 0.100 

Abell 548b^ 0.042 

Abell 746 0.232 

Abell 1240 0.159 

Abell 1612 0.179 

Abell 1664 0.128 

Abell 2061 0.078 

Abell 2163 0.203 

Abell 2255 0.081 

Abell 2256 0.059 

Abell 2345 0.176 

Abell 2744 0.308 

Abell 3365 0.093 

Abell 3376 0.046 

Abell 3667 0.055 

Coma cluster 0.023 

CIZA J2252.8+5301 0.192 

MACSJ0717.5+3745 0.555 

RXC J1053.7+5452 0.070 

RXCJ1314.4-2515 0.244 

ZwCl 0008.8+5215 0.104 

ZwCl 2341.1+0000 0.270 



45 

6 

23 

87(A), 40(B) 

25 

13(N), 15(S) 

62 

3 

15 

1 

61 

26 

13(W), 24(E) 

20 

20(W), 141(E) 

3(W), 8(E) 

13(W), 0(E) 

22 

27 (N), 4(S) 

2 

19 

30(W), 99(E) 

10(W), 11(E) 

14(N), 8(S) 



78 

64 

39 

84(A), 48(B) 

56 

85(N), 85(S) 

69 

78 

77 

67 

38 

77 

87(W), 43(E) 

81 

82(W), 43(E) 

85(W), 67(E) 

75(W), 75(E) 

60 

67(N), 84(S) 

79 

76 

26(W), 46(E) 

76(W), 85(E) 

88(N), 63(S) 



for double relics two values are measured (N=north, S=south, E=east, W=west) 

references: Abell 225 5 (main relic only) |Pizzo &de~B ruyn (2009]>; [Govoni et al.|p005]>, MACS J07 17.5+3 745 |van Weeren eTaIl(l2009c| ); 
Bonafede et al](|2009a|), Abell 548b (2 relics) |Feretti et al.|(|2006| 7Abell 2256 Clarke & Ensslin (2006 ), Abel l 521|Giacintucci et a l. (2008 ), Coma 



Cluster Giovannini et al 



1991 



Abell 1664 Govoni et al. 



2001 



Abell 3376 Bagchi et al. 



2006 



; |Brown & Rudnick[( [201l] )7Abell 2163 |Feretti et al.| ( |2001| ), Abell 2744 |Govoni et al.|([20011), Abell 115 and 
h. Symmetric double relics incl uded from the literature are: Abell 1240 and Abell 234 5 Bonafe de" et al.| ([2009b), 
_J|, Abell 3 667|Rottgering et al.|(T997]>, ZwCl 0008.8+5215 [van" Weeren et al.| ( |2011a|), CIZA J2242.8 +5301|van 
IWeeren et al.| ( |2010| >, ZwCl 2341. l+0000 |van Weeren et al.|p009a) ; |Giovannini et al.| ( |2010| ), RXC Jl 3 14.4-25 15 |Venturi et al| ( [2007l ). 

a identification uncertain (either relic or halo) 

^A548b: A and B refer to the names used in |Feretti etall ( |2006| ) 



used in the simulation the radio emission strongly increases with 
Mach number (Hoeft & Briiggen 2007 ), a rather small increase 
of the Mach numbers would have a large effect on the result- 
ing radio luminosity. (We note that the simulation of |Hoeft & 
Briiggen (2007 ) is also based on DSA.) Finally, the merger rate 
increases with redshift. However, it needs to be clarified from 
simulations why an increase of the radio power with X-ray lu- 
minosity is expected. The fractions of relics in X-ray selected 
cluster samples is therefore a powerful tool to put constraints on 
the models used in the simulations, and hence on the evolution 
of magnetic fields and on particle acceleration in the ICM. 



4.2. X-ray peak and galaxy distribution separation 

Cluster mergers are thought to decouple the baryonic matter 
component from the dark matter (DM). This causes and offset 
between the gravitational center (measured from lensing) and X- 
ray center of the cluster. A clear example is the "Bullet cluster" 
(e.g., |Clowe et a l. 2006 ), but the effect has also been measured 
for other clusters (e.g., Shan et al.|2010| ). Radio relics and halos 
are mostly found in merger cluster and are thus good candidates 
to measure this effect. We do not have lensing measurements for 
the clusters in our sample, but simply taking the galaxy distri- 
bution as a proxy for the dark matter distribution we note off- 
sets between the X-ray emission from ROSAT and the galaxy 
distribution in the clusters A2034, A3365, A2061, A523 and 




Fig. 23. Histograms showing the distribution of the major and 
minor axis ratio. The red histogram is for clusters from the 
NORAS and REFLEX surveys, the black histogram is for clus- 
ters containing giant radio relics (see the caption of Fig. [22]). The 
black histogram was scaled by a factor of six for easier compar- 
ison with the NORAS-REFLEX sample. Clusters for which the 
2 dimensional Gaussian fit did not converge were not included. 
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Table 5. Clusters with radio relics in the NORAS and REFLEX surveys 



Cluster name 



^i.4ghz (total for cluster) L x , 0.1-2.4 b 
lO^WHz" 1 lO^ergs- 



A3376 


0.0468 


1.6 


A2744 


0.3066 


7.7 


A4038 c 


0.0292 


0.1 


A0548W^ 


0.0424 


0.5 


RXCJ1314.4-2515 


0.2439 


8.2 


A133 c 


0.0569 


1.1 


A1300 


0.3075 


6.3 


A13 c 


0.0940 


0.9 


A2345 


0.1760 


6.2 


A521 


0.2475 


3.4 


A754 


0.0542 


0.04 


A85 c 


0.0555 


0.3 


A2163 


0.2030 


2.8 


A1612 


0.1797 


7.9 


A523* 


0.1034 


1.8 


Coma cluster 


0.0231 


0.29 


A781^ 


0.2952 


5.8 


A2034 a 


0.1130 


1.17 


A1758 


0.2799 


4.1 


A746 d 


0.232 


6.8 


RXCJ1053.7+5452 


0.0704 


0.2 


A2255 


0.0809 


0.8 


A2256 


0.0581 


4.2 



1.01 

11.68 

1.00 

0.10 

9.92 

1.40 

12.12 

1.24 

3.91 

7.44 

3.79 

5.18 

19.62 

2.41 

0.89 

3.63 

4.15 

3.56 

10.90 

3.68 

0.44 

3.08 

3.69 



references: Abell 754 Macario et al. (2011); Abell 1758 |Giovannini et al.l([2009] ); Abell 781 |Venturi et"aH ( [20TT] ); Abell 1300 |Reid et aTR l999); 

Abell 13, 85, 133, and 4038 Sl ee et al.|(2Ml) ; for other clusters see Table|4] 

a identification uncertain (either relic or halo) 

b below X-ray flux limit of F x , 0.1-2 a kev = 3 X 10~ 12 erg s _1 cm -2 

c radio phoenix, not included in sample 




0.2 0.3 

redshift 



Fig. 24. Lx-redshift distribution for the NORAS and REFLEX 

surveys. The solid blue line is the flux cutoff of 3.0 x 
10" 12 erg s" 1 cm -2 we use for comparison with the relic clus- 



ter sample. 



RXC J1053. 7+5452 for example. These clusters would there- 
fore make interesting targets for future X-ray and lensing obser- 
vations. 
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Fig. 25. P 1.4GHz -redshift distribution for the relics in Table[5j The 
solid line is for a relic flux limit of 10 mJy, the dashed line for 
30 mJy. Radio phoenices are not included. 



5. Conclusions 

We have presented WSRT, VLA and GMRT observations of 
galaxy clusters with diffuse emission selected from the WENSS 
and NVSS surveys. We find peripheral radio relics in the clus- 
ters Abell 1612, Abell 746 and CIZA J0649.3+1801 and a 
smaller relic Abell 2034. Abell 3365 seems to host a double ra- 
dio relic system. Our observations reveal radio halos in the clus- 
ters Abell 746 and RX J0107. 8+5408. We confirm the presence 
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r 15 AA 2 ° n ™ w 0.0 0.1 0.2 0.3 

L x [I0 44 erg/sj redshift 

Fig. 26. Histograms showing the X-ray luminosity (left) and redshift (right) distribution. Blue histograms shows the NORAS- 
REFLEX sa mple, red the relic cluster sample. The solid blue line displays the predicted luminosity/redshift distributions from Nuza 
et al. ( |201l| for clusters with fluxes > 3.0 x 10" 12 erg s" 1 cm -2 , while the solid red line is the prediction for clusters hosting relics 
in the simulation. The fraction of clusters with relics is given by the black solid line (the ratio of the blue and red lines). 



of radio relics in Abell 2061, RXC J1053.7+5452 and diffuse 
emission in Abell 523 and Abell 2034 (for which the classifica- 
tion is uncertain). In addition, we detect the radio halo Abell 697 
providing additional flux measurements around 1.4 GHz. 

By constructing a sample of 35 relics we have found that 
relics are generally located along the major axis (which can 
be used as a proxy for the merger axis) of the cluster's elon- 
gated ICM. Their orientation is mostly perpendicular to this ma- 
jor axis. The distribution of the major-minor axis ratios for the 
relic cluster sample is broader than that of the NORAS -REFLEX 
sample. These results are consistent with the scenario that relics 
trace shock waves which form along the merger axis of clus- 
ters. We also compared the redshift and X-ray luminosity dis- 
tributions of clusters with relics to a sample of clusters from 
the NORAS and REFLEX surveys. Interestingly, we find indi- 
cations that the observed fraction of clusters hosting relics in- 
creases with X-ray luminosity and redshift. However, selection 
biases for radio relics play a role and this needs to be investi- 
gated further. The significantly improved sensitivity of upcom- 
ing radio telescopes (e.g., the EVLA and LOFAR) will allow to 
find many more radio relics. Correlating this larger relic sample 
with cluster samples selected at X-rays wavelengths will pro- 
vide a powerful tool to constrain the evolution of relics in galaxy 
clusters. 
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